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Abstract 

Statements concerning Extended Unitarity (EU) are clarified. 
PACS numbers: 11.80.Et, 1360.6e, 14.40. Cs. 

The objective of this Comment is to expose points of uncertainty and disagreement concerning 
Extended Unitarity (EU). Liu et al. PQ make several assertions criticising an earlier publication 
of mine [2], where fits to four sets of data are compared with the form of Extended Unitarity 
formulated by Au, Morgan and Pennington [3]. The principles should apply to the overlap of 
any two resonances having the same quantum numbers. Both my publication and that of Liu et 
al. concern a and /o(980) where they overlap near 1 GeV. Here a is a shorthand for the broad 
7i7i component going through a phase shift of 90° in the mass range 900-940 MeV. 

For 7i7i elastic scattering, a Breit-Wigner resonance has an elastic scattering amplitude 
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p is Tin phase space and N(s) is real. Watson's theorem [1] amounts to the fact that the Breit- 
Wigner denominator D(s) and the associated pole are universal. The numerator is not universal. 

Fig. 1(a) shows the tttt mass projection from BES II data for J/ty — > ujtt + 7t~ [5]. The peak 
near 500 MeV may be fitted with the a pole taking N(s) = 1, unlike elastic scattering where 
N(s) is constrained by unitarity to be Mr„(s). There may be a small dip in the data due to 
/o(980) in the bin just above the KK threshold. The a amplitude has a different s-dependence 
to elastic scattering and the relative magnitude of /o(980) to a is also different. 

Let us denote the tttt production amplitude by Ak(s), where k refers to channels tttt, KK, 
etc. Liu et al. use the relation 

Im(A0 = -p k T* k A k , (3) 

where T ik is the 2 — > 2 amplitude. The form of EU used by Liu et al. is that the phase of T is 
identical to that of elastic scattering. They do not specify if this is the phase shift 5 or the phase 
(j) measured from the bottom of the unitarity circle. These differ above the inelastic threshold. 
The Omnes relation [6] connects the magnitude of the elastic scattering amplitude to <ft using 
analyticity, so <fi appears to be the most relevant angle and will be used here. 

The critical assumption of EU is that the relative phases of o and /o in T must be identical 
in production reactions and elastic scattering. This is an assumption going beyond Watson's 
theorem [I]. 

In elastic scattering, a and fo amplitudes are constrained below the KK threshold to move 
around the unitary circle on the Argand diagram. In Cern- Munich data [7] , it is obvious by eye 
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Figure 1: (a) The tttt mass projection for BES II data on J/$ ^ um + 7i^; the upper histogram 
shows the fit to EU in the form of AMP, and the lower histogram the irn S-wave intensity; (b) 
an enlarged view of the corner of the Dalitz plot; cross-hairs show the point of intersection of 
/o(980) with the strong &!(i235) component. 



that phase shifts of a and /o(980) add to a good approximation near 1 GeV. This additivity of 
phases may be accomodated by multiplying 5"- matrices r]exp(2i5), though it is an open question 
whether rj parameters should multiply. For present considerations concerning the overlap of a 
and fo, this is not critical, since the inelasticity of the a rises fairly slowly at the KK threshold 
because the a has a large Tin width. 

For a production process the relative magnitudes of a and f are different, because of different 
matrix elements connecting initial and final states. It is not obvious that T of different resonances 
must combine in the same way as for elastic scattering. In a 2 — > 2 process, a and fo amplitudes 
combine to make an overall amplitude T which obeys the relation 

ImT = TT*. (4) 

Inal— >3or2— >3 process, the boundary condition for the initial state is different. Why 
can /o(980) not be completely absent in some production processes? In that case, Eq (4) is still 
valid for the a amplitude alone, but EU demands that the phase of /o(980) still appears in the 
production process despite its magnitude being absent. Note that the process J/\& — > U7i + n^ 
accounts for only 1% of all J/\l/ decays. 

These questions were much discussed in the early 1960s and led to the isobar model, where 
each pole is multiplied by complex coupling constants and amplitudes are added, not multiplied. 
For experimentalists, this is a simple form with which to parametrise data, except possibly in the 
small region where they overlap and EU could play a role. Even there it parametrises average 
phase differences between resonances. In the production process, the matrix element involves 
an unknown integral between the initial and final-state wave functions. This is what Liu et al. 
seek to model. There is no dispute that multiple scattering processes do affect phases. 

The Omnes relation, based on analyticity, allows the magnitude of the tttt amplitude to be 
derived from phases (p. A subtlety is that the magnitude can be multiplied by a polynomial. This 
polynomial accomodates, for example, a form factor arising from the overlap integral between 
initial and final states. Liu et al. note this polynomial, but do not give details of whether or 
how they use it. 
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On Fig. 1(a). the magnitude of the a amplitude is not the same as for elastic scattering. 
My approach, consistent with the data within errors, is to take N(s) = 1. Liu et al. may use 
instead an explicit polynomial relating elastic scattering and production, but they do not say. 
They derive relative magnitudes of a and fo from equations given in earlier work by Lahde 
and Meissner [8]. This takes a and fo to be nn and ss linear combinations, though results 
would be similar assuming 4-quark compositions (apart from SU3 Clebsch-Gordan coefficients). 
The matter of different SU3 Clebsch-Gordan coefficients is not presently a critical issue for EU, 
though it would bear upon the relation of a and fo to 2-quark, 4-quark and/or meson- meson 
components in wave functions. If the relative magnitudes of /o(980) and a are the same in 
J/^f — > WK7T data and elastic scattering, the polynomial is 1 (except for standard form factors), 
but if the relative magnitudes are different, a polynomial is required peaking at the mass of 
/o(980) and reproducing its shape. What does this polynomial mean? 

From the point of view of the isobar model, this is inconvenient. The idea of the isobar model 
is that form factors should have small and slowly varying effects, usually negligible. Liu et al 
do not show the Argand diagram for their amplitude. It needs to be exposed what polynomial 
P(s) is being used, so the reader knows what price is being paid to preserve EU. 

Paradoxes can arise in limiting cases. Suppose the f amplitude is strictly zero in a production 
process. One then sees the a amplitude alone. According to Watson's theorem, the production 
amplitude must then have the phase of the a. The hypothesis of EU requires a polynomial which 
exactly cancels the magnitude and phase of the fo in elastic scattering; but the cancellation of 
the phase is not allowed, therefore a zero fo amplitude appears to be forbidden by EU. My 
opinion is that this is a critical defect in the hypothesis of EU. 

Let us return to details of J/^> —> unit data. There were two fits reported in the BES 
publication [5]. The second does contain a large /o(980) signal. However, this gives a fit worse 
in x 2 by over 200 compared to a free fit to / (980). For fine details, the reader should consult 
Ref. [2], but Fig. 1 presents essential points. Fig. 1(b) shows an enlargement of the corner of 
the Dalitz plot where the f crosses the strong &x(1235) band in the data; the intersection is 
shown by the cross-hairs. If there were a strong fo component with the same phase as a, one 
should see the interference dip between them as a diagonal band running at 45° from top left to 
bottom right. It should show a distinctive interference as it crosses the 6i (1235) band centred 
on the vertical cross-hair. This interference should have a full width at half-maximum of 0.07 
GeV 2 in m 2 (co>7r). There is no sign of this interference with &!(1235), o or / 2 (1270). This is why 
the fitted /o(980) components is very small in my fit to the data, Ref. [2]. 

Liu et al. do not have access to the Dalitz plot. Their Fig. 5 shows structure just below 1 
GeV which seems to require an fo at least as large as the a amlitude in this mass region. Firstly 
it is not obvious how to reconcile this with J/\I> — > 07r + 7r~ data where the a amplitude is much 
smaller than the strong / (980) peak. Secondly the precipitious drop in the amplitude fitted to 
ujix+ix~ data just below the KK threshold should generate a cliff-like diagonal structure on Fig. 
1(b). The fit of Liu et al. is inconsistent with the isobar model fit shown as the dashed curve 
on Fig. 1(a). The broad diagonal band of Fig. 1(b) arises from interferences of a with &i (1235) 
and / 2 (1270). 

It would help the discussion if the BES II collaboration would make publicly available the 
data set produced by myself and collaborators, the corresponding Monte Carlo data and the 
sidebin events used to subtract the 14% experimental background. The decay plane of the u 
provides delicate information defining accurately the &i(1235)7r signal, which then serves as an 
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Figure 2: (a) AFS data compared with the isobar model fit; (b) enlargement near 1 GeV: the full 
curve is the isobar model fit, the chain curve shows what appears in elastic tttt scattering, and 
the dashed curve shows the isobar model fit supplemented by a freely fitted I = 2 tttt amplitude. 



interferometer determining magnitudes and phases of amplitudes crossing it. 

A much clearer test of EU arises in data from the ISR [9j for central production of tttt in 
the process pp — > pp(tttt). The central tttt pair is far removed in rapidity from final states 
protons, and there is no evidence for a Deck effect, i.e. resonant pir or pirn combinations in the 
central region. In this case, EU predicts that relative phases of o and fo should be identical to 
2^2 processes. The data are easily fitted by the isobar model. It requires relative intensities 
/o(980)/cr ~ 60% of elastic scattering at the KK threshold. The fitted /„ phase is (57 ± 7)° 
below the EU prediction. Features of the fit are displayed in Figs. 9 and 10 of Ref. [2]. Fig. 2 
summarises essential points. 

Fig. 2(a) shows the isobar model fit (which includes a small slowly varying form factor to 
reproduce the precise shape of the sigma pole up to 0.8 GeV). Fig. 2(b) shows details of the 
mass range near 1 GeV. The full curve is the isobar model fit without any I = 2 tttt amplitude. 
The chain curve shows what is observed in elastic tttt scattering, suitably normalised. These 
two curves are very different in the region of /o(980), hence the phase difference required in the 
isobar model fit. There could be an / = 2 production amplitude. The dashed curve shows the 
result of fitting this freely. The fit is still poor and the 1 = 2 amplitude is twice as large as in 
elastic scattering. My view is that these data rule out EU in the form of Liu et al. 

Morgan and Pennington [TU] did manage to fit central production data with EU but required 
an additional third-sheet pole at M = 978— i28 MeV. Since then, data on /o(980) have improved, 
and there is no sign of this additional narrow third-sheet pole in BES II data. 

In Ref. [2], an empirical relation is presented which does fit the ISR data. The assumption 
is made that the production process can be considered as a black box. Then it is postulated 
that whatever linear combination of a and fo is produced will rescatter according to the usual 
relation Im T = TT*, i.e. the pions rescatter as an isolated pair with magnitudes different from 
elastic scattering. If the asymptotic combination is written T = T a + f3Tf exp 2 ** within the 
elastic regime, and f3 is real, it can be shown that 

2^ = 2<5 ff -<5 /o +sin- 1 (/3sin ( 5 /o ). (5) 

If /? ^ 1, this is a different relation from the 2 — > 2 process. As (3 decreases from 1, \I/ falls 
rapidly. In Ref. [2], the exact form Im T = TT* is used, including effects of coupling to KK. 
Although Liu et al. dismiss this approach, it does have the merit of correctly predicting the 
observed phase difference between a and fo within errors. 
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No analyses so far include the dispersive cusp at the KK threshold in both a and /o(980). 
The cusp originates from a discontinuity in slope of the KK amplitude at threshold. This 
causes a corresponding discontinuity in slope of the real part of the amplitude, hence the cusp 
[TT] . I have made some preliminary fits to BES II data including this cusp, and can offer some 
cautionary comments to others who may try something similar. In order to make the dispersion 
integral of Eq. (2) converge, a form factor must be included in the KK channel. Whatever 
form factor is adopted, there is a narrow peak in the real part of the hit amplitude exactly at 
the KK threshold. The narrow width implies that mass resolution is of critical importance in 
any comparison with data. The width of the cusp is quite stable. However, its height above 
zero real part does depend on exactly what form factor is used. Either this form factor must 
be predicted from a model of /o(980) or the height of the cusp must be treated as a variable in 
fitting data. This source of uncertainty adds to difficulties in deciding the nature of /o(980). 

The advice to experimentalists is to use the isobar model, which is simple and well defined. 
It takes account empirically of phases such as Liu et al. are trying to fit in a different way. If 
mass resolution is very good, it will be necessary to include some s-dependence in T n7T similar 
to the cusp. Theorists can then compare the fit with models of f and the cusp those models 
predict. It may well be necessary to allow for f —> r\r\. 

In conclusion, the free fit made by Liu et al. to intermediate steps in J/ty — > oj-k+h^ 
accounts for phases used in the isobar model, but the magnitude of /o(980) signal they fit 
looks inconsistent with Fig. 1(b). The ISR data are in strong disagreement with EU. In those 
data, there should be no significant perturbations due to other resonances and the production 
mechanism is clean. 

References 

B. Liu et al., arXiv: 0901.1185 
[2] D.V. Bugg, Eur. Phys. J C54 (2008) 73 [arXiv: 0801.1908] 
[3] K.L. Au, D. Morgan and M.R. Pennington, Phys. Rev. D35 (1987) 1633 
[4] K.M. Watson, Phys. Rev. 88 (1952) 1173 

[5] M. Ablikim et al. (BES II Collaboration) Phys. Lett. B598 (2004) 149 [arXiv: 



[6] 
[7] 
[8 
[9 
[10 

[11 



hep-ex/0406038] 



R. Omnes, Nu. Cim. 8 (1958) 316 

B. Hyams et al, Nucl. Phys. B64 (1973) 134 

T.A. Lahde and U.-G. Meissner, Phys. Rev. D74 (2006) 034021 [arXiv: |hep-ph/0606133] 

T. Akesson et al., Nucl. Phys. B264 (1986) 154 

D. Morgan and M.R. Pennington, Phys. Rev. D48 (1993) 11851 

D.V. Bugg, J. Phys. G: Nucl. Phys. 35 (2008) 075005 [arXiv: 0802.0934] 



5 



